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We examine the effect of interactions between the eiectrons on the Landauer-Biittiker conductances 
of some systems of quantum wires with different geometries. The systems include a long wire with 
a stub in the middle, a long wire containing a ring which can enclose a magnetic flux, and a system 
of four long wires which are connected in the middle through a fifth wire. Each of the wires is taken 
to be a weakly interacting Tomonaga-Luttinger liquid, and scattering matrices are introduced at all 
the junctions present in the systems. Using a renormalization group method developed recently for 
studying the flow of scattering matrices for interacting systems in one dimension, we compute the 
conductances of these systems as functions of the temperature and the wire lengths. We present 
results for all three regimes of interest, namely, high, intermediate and low temperature. These 
correspond respectively to the thermal coherence length being smaller than, comparable to and 
larger than the smallest wire length in the different systems, i.e., the length of the stub or each 
arm of the ring or the fifth wire. The renormalization group procedure and the formulae used to 
compute the conductances are different in the three regimes. In particular, the dimensionality of the 
scattering matrix effectively changes when the thermal length becomes larger than the smallest wire 
length. We also present a phenomenologically motivated formalism for studying the conductances 
in the intermediate regime where there is only partial coherence. At low temperatures, we study 
the line shapes of the conductances versus the energy of the electrons near some of the resonances; 
the widths of the resonances are found to go to zero with decreasing temperature. Our results show 
that the Landauer-Biittiker conductances of various systems of experimental interest depend on the 
temperature and lengths in a non-trivial way when interactions are taken into account. 
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I. INTRODUCTION 

The increasing sophistication in the fabrication of 
semiconductor hctcrostructures and carbon nanotubes in 
recent years have made it possible to study electronic 
transport in different geometries. For instance, three- 
arm and four-arm quantum wire systems have been fabri- 
cated by voltage-gate patterning on the two-dimensional 
electron gas in GaAs heterojunctions [1,2]. Other sys- 
tems of interest include Y-branched carbon nanotubes 
[3], crossed carbon nanotubes [4], mesoscopic rings [5,6], 
and quantum wire systems with stubs [7]. There have 
also been many theoretical studies of transport in sys- 
tems with various geometries [8-12]. 

Studies of ballistic transport in a quantum wire (QW) 
have led to a clear understanding of the important role 
played by both scattering of the electrons and the in- 
teractions between the electrons inside the QW [13-16]. 
The scattering can occur either due to impurities inside 
the QW or at the contacts lying between the QW and its 
reservoirs. A theoretical analysis using bosonization [17] 
and the renormalization group (RG) method typically 
shows that repulsive interactions between electrons tend 
to increase the effective strength of the back-scattering 
as one goes to longer length scales; experimentally, this 



leads to a power-law decrease in the conductance as the 
temperature is reduced or the wire length is increased 
[18]. Motivated by this understanding of the effects of 
interaction on scattering, there have been several studies 
of the interplay between the effects of interactions on one 
hand, and either a single junction between three of more 
QWs [19-22], or more complicated geometries [23-26] on 
the other. Using a RG technique introduced in Ref. [27], 
the effects of a junction (which is characterized by an 
arbitrary scattering matrix S) has been studied in some 
detail [21]. It is now natural to extend these studies to 
systems of QWs which are of experimental interest and 
which can have with more complicated geometries involv- 
ing more than one junction. 

In this paper, we will study the effect of interactions on 
the Landauer-Biittiker conductances of three systems of 
quantum wires with different geometries. These systems 
are shown in Figs. 1-3, and we will refer to them as 
the stub, the ring and the four-wire system respectively. 
The stub system consists of two long wires, labeled as 
1 and 3, with a stub labeled as 2 being attached to the 
junction of 1 and 3. The ring consists of two long wires, 
labeled as 1 and 3, between which there is a ring which 
can possibly enclose a magnetic flux; the two arms of the 
ring, labeled as 2 and 4, will be assumed to have the same 
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length for convenience. The four-wire system consists of 
four long wires labeled as 1, 2, 3 and 4. The junction 
of 1 and 2 is connected to the junction of 3 and 4 by a 
fifth wire labeled as 5. The length of wire 2 in the stub 
system, the length of each of the arms 2 and 4 in the 
ring system, and the length of wire 5 in the four-wire 
system will all be denoted by L$- Each of the junctions 
present in the different systems is governed by a 3 x 3 
scattering matrix S which is unitary. We will assume that 
each of the wires in the various systems can be described 
as a one-channel weakly interacting Tomonaga-Luttinger 
liquid (TLL). For simplicity, we will ignore the spin of 
the electrons in this paper. 
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FIG. 1. The stub system, showing two long wires labeled 
as 1 and 3, and a stub labeled as 2. The lower end of the 
stub where three wires meet and the upper end of the stub 
are denoted by A and B respectively. 
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FIG. 2. The ring system, showing two long wires labeled 
as 1 and 3, the two arms of the ring labeled as 2 and 4, and 
two three-wire junctions labeled as A and B. 




FIG. 3. The four-wire system, showing four long wires la- 
beled as 1, 2, 3 and 4, a connecting wire in the middle labeled 
as 5, and two three-wire junctions labeled as A and B. 



of the S'-matrices at the various junctions of the differ- 
ent systems. In Sec. Ill, we will describe the procedure 
for computing the transmission probabilities (and con- 
ductances) of a system given the form of the S'-matrices 
at all its junctions. It turns out that both the RG pro- 
cedure and the route from the S-matrix to the conduc- 
tances depend on the range of temperatures that one is 
considering. There is a length scale, called the thermal 
coherence length Lt, which governs the typical distance 
beyond which the phase of the electron wave function be- 
comes uncorrelated with its initial phase. The regimes of 
high, intermediate and low temperatures are governed re- 
spectively by the condition that Lt is much smaller than, 
comparable to or much larger than the length scale L$ 
defined above for the three systems; correspondingly, we 
have complete incoherence, partial coherence and com- 
plete coherence for the phase. The intermediate tem- 
perature range is the most difficult one to study, both 
for using the RG method and for computing the con- 
ductances. Based on some earlier ideas [28,29], we will 
describe a phenomenological way of introducing partial 
coherence which will lead to expressions for the trans- 
mission probabilities which interpolate smoothly between 
the coherent and incoherent expressions. 

In Sees. IV- VI, we will apply the formalism outlined 
in the previous sections to the stub, ring and four-wire 
systems respectively. In each case, the transmission prob- 
abilities at intermediate and low temperatures (i.e., the 
partially and completely coherent regimes) will be found 
to depend sensitively on the phase ij — e l2kFLs ; here Uf 
is the wave number of the electrons which are assumed 
to come into or leave the QW system with a momen- 
tum which is equal to the Fermi momentum in the reser- 
voirs. In particular, certain values of rj can lead to reso- 
nances and anti-resonances i.e., maxima and minima in 
the transmission probabilities. In the ring system, there 
is another important phase which governs the possibility 
of resonance, namely, e le( t>B/fic^ ^grg ^ B j s ^he magnetic 
flux enclosed by the ring, and e and c are the electron 
charge and the speed of light respectively. In each sys- 
tem, we will see how the conductances vary with the 
temperature in a non-trivial way as a result of the inter- 
actions. This is the main point of our paper, namely, that 
interactions between the electrons lead to certain power- 
laws in the temperature and length dependences of the 
conductances of experimentally realizable quantum wire 
systems. 



II. RENORMALIZATION GROUP METHOD FOR 
SYSTEMS WITH JUNCTIONS 



In Sec. II, we will first summarize the RG method 
developed in Ref. [21] for studying the flow of the S- 
matrix at a junction due to the interactions in the dif- 
ferent wires connected to that junction. We will then 
describe our method for carrying out the RG analysis 



In this section, we will first present the RG procedure 
developed in Ref. [21] for studying how the effect of a 
single junction varies with the length scale. We will then 
describe how the RG method has to be modified when a 
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system has more than one junction. 

A junction is a point where N semi-infinite wires meet. 
Let us denote the various wires by a label i, where 
i = 1, 2, • • • , N. As we approach the junction, the incom- 
ing and outgoing one-electron wave functions on wire i 
approach values which are denoted by ipn and ipoi re- 
spectively; we can write these more simply as two N- 
dimensional columns ipi and ipo- The outgoing wave 
functions are related to the incoming ones by a N x N 
scattering matrix, 



ipo = S ipi 



(1) 



Current conservation at the junction implies that S must 
be unitary. (If we want the junction to be invariant under 
time reversal, S must also be symmetric). The diagonal 
entries of S are the reflection amplitudes ra, while the 
off-diagonal entries are the transmission amplitudes tij 
to go from wire j to wire i. We will assume that the 
entries of S do not have any strong dependence on the 
energy of the electrons. 

Let us now consider a simple model for interactions 
between the electrons, namely, a short-range density- 
density interaction of the form 



Hir,* — — 



dxdy p(x) V(x - y) p(y) , (2) 



where V(x) is a real function of x, and the density p{x) is 
given in terms of the second-quantized fermion field ^>(x) 
as p = <J't\I>. [The assumption of a short-ranged inter- 
action is often made in the context of the TLL descrip- 
tion of systems of interacting fermions in one dimension.] 
We define a parameter g 2 which is related to the Fourier 
transform of V(x) as 52 = V(0) — V(2kp)- Different wires 
may have different values of this parameter which we will 
denote by gn. For later use, we define the dimensionless 
constants 
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where we assume that the velocity 



vf = 



hkr. 



m 



(3) 



(4) 



is the same on all wires. In this work, we will be in- 
terested in the case in which the interactions are weak 
and repulsive, i.e., the parameters aj are all positive and 
small. 

We are now ready to present the RG equation for the 
matrix S which was derived in Ref. [21]. It will be useful 
to briefly discuss the derivation of the RG equation. A 
reflection from a junction, denoted by the amplitude ra in 
wire i, leads to Friedel oscillations in the electron density 
in that wire. If x denotes the distance of a point from the 
junction, the form of the oscillation at that point is given 
by the imaginary part of rue FX / '(2nx). As a result of 



the interactions, an electron traveling in that wire gets 
reflected from these oscillations. The amplitude of the 
reflection from the oscillations is proportional to anru/2 
if the electron is reflected away from the junction, and to 
a^r^/2 if the electron is reflected towards the junction. 
These reflections renormalize the bare S'-matrix which 
characterizes the junction at the microscopic length scale. 
The entries of S therefore become functions of the length 
scale L; we define the logarithm of the length scale as 
I = \n(L/d), where d is a short-distance cutoff such as the 
average interparticle spacing. It is convenient to define a 
N x N diagonal matrix M whose entries are given by 
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Then the RG equation for S is found to be [21] 

^ = M - SM^S 
dl 



(5) 



(6) 



to first order in the a*. (This equation is therefore per- 
turbative in the interaction strength). One can verify 
from Eq. (6) that S remains unitary under the RG flow; 
it also remains symmetric if it begins with a symmet- 
ric form. The fixed points of Eq. (6) are given by the 
condition S'M t = MS\ i.e., SM* must be Hermitian. 

We can study the linear stability of a fixed point by 
deviating slightly from it, and seeing how the deviation 
grows to first order under the RG flow. Let us denote 
a fixed point by the matrix So and a deviation by eS±, 
where e is a small real parameter and Si is a matrix; we 
require that the matrix S = Sq + eSi is unitary up to 
order e. (We can think of Si as defining the "direction" 
of the deviation). We substitute S in Eq. (6) and then 
demand that Si should take such a form that the RG 
equation reduces to 



dl 



pe 



(7) 



where p is a real number. We then call the direction Si 
stable, unstable and marginal (to first order) if p < 0, > 
and respectively. All fixed points have at least one ex- 
actly marginal direction which corresponds to multiply- 
ing the matrix So by a phase; clearly this leaves Eq. (6) 
invariant. 

In this paper, we will be concerned with the RG flow 
of S- matrices which are 2, 3 and 4 dimensional. For 
convenience, we will assume certain symmetries in each 
of these cases. It is useful to discuss these symmetries 
here, and how they lead to some simplifications for the 
RG flows. 

We first consider a two-wire system in which there is 
complete symmetry between the wires which we will label 
as 1 and 3. Namely, the interaction parameters are equal, 
011 = «3 = a, and the scattering matrix has the form 



S2D — 



a b 
b a 



(8) 
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Unitarity implies that we can parametrize a and b as 



iXe" 



b = 



1 + iX ' 



1 + iX ' 



(9) 



where A and 9 are real. Eq. (6) then leads to the following 
differential equations 



dX 

~dl 
d9_ 



= a X , 
aA 



1 + A 2 



(10) 



The reflection and transmission probabilities \a\ 2 and |6| 2 
only depend on A. For a > 0, we see that there is an 
unstable fixed point at A = 0, and a stable fixed point at 
A = oo. If A is not zero initially (i.e., at the microscopic 
length scale d), then it flows to infinity at long distances. 
Hence b goes to zero as 



t 



(11) 



a approaches 1, and the two wires effectively get cut off 
from each other. This is in agreement with the results 
obtained using bosonization [17]. 

Next, we will consider the 3x3 case. Here we will 
assume that there is complete symmetry between two of 
the wires, say, 1 and 2, and that the 5-matrix is real. 
Namely, ct\ — ct2, and S takes the form 



S3D — 




(12) 



where r', t' and t are real parameters which, by unitarity, 
satisfy 



1 + r' , 
- 1 - 2/ , 



t=V(-2r>) (1 + r') , 



(13) 



and — 1 < r' < 0. The RG equations in Eq. (6) can be 
written purely in terms of the parameter r 1 as 



r' (1 + r 1 ) [air 7 + a 3 (l + 2r / ) 



(14) 



If ai, 0:3 > 0, we have stable fixed points at r' = (where 
there is perfect transmission between wires 1 and 2, and 
wire 3 is cut off from the other two wires) and — 1 (where 
all three wires are cut off from each other) . There is also 
an unstable fixed point at 



a 3 



(15) 



If r' starts with a value which is greater than (or less 
than) this, then it flows to the value (or —1) at large 



distances. [We should point out that the fixed point in 
which wires 1 and 2 transmit perfectly into each other 
and wire 3 is cut off is stable only within the restricted 
space described by Eqs. (12-13). If we take a general 
unitary matrix S 3 d, then this is not a completely stable 
fixed point. The only stable fixed point in the general 
case is the one in which all three wires are cut off from 
each other [21].] 

Finally, let us consider the 4x4 case. Here we will 
be interested in a situation in which there is complete 
symmetry between wires 1 and 2, and and between wires 
3 and 4; further, we will take the values of a, in all the 
wires to be equal to a. The 5-matrix takes the form 



Sad = 



( a b c c\ 

b a c c 

c c a b 

\c c b a J 



(16) 



where a, b and c are all complex. Unitarity implies that 
these parameters can be written in terms of three in- 
dependent real variables. There does not seem to be 
a convenient parametrization in terms of which the RG 
equations take a simple form. We therefore have to study 
the RG equations in Eq. (6) numerically; the results will 
be described in Sec. VI. However, the fixed points of the 
RG equations and their linear stabilities can be found 
analytically. There are three kinds of fixed points. 

(i) \a\ = 1, and b = c = 0. This corresponds to all the 
wires being cut off from each other. This fixed point is 
stable in two directions, and is exactly marginal in one 
direction (corresponding to a phase rotation of a). 

(ii) |&| = 1, and a = c = 0. This corresponds to perfect 
transmission between wires 1 and 2, and between wires 
3 and 4, but no transmission between any other pair of 
wires. This fixed point is unstable in one direction (where 
it flows to the fixed point described in (i)), and marginal 
in two directions. One of these marginal directions turns 
out to be unstable at a higher order, and the RG flow 
eventually takes it to the third fixed point described be- 
low. The other marginal direction corresponds to a phase 
rotation of b. 

(iii) |a| = 1/2, b = —a, and c = ±a. This is a spe- 
cial point which corresponds to the maximum possible 
transmission with complete symmetry between all the 
four wires. This fixed point is unstable in one direction 
(where it flows to the fixed point in (i)), stable in a sec- 
ond direction (where it flows in from the fixed point in 
(ii)), and exactly marginal in the third direction (corre- 
sponding to a simultaneous phase rotation of a, b and c). 
The fact that (iii) is stable in one direction and unsta- 
ble in another, means that an interesting cross-over can 
occur as a result of the RG flow. Namely, one can begin 
near (ii), approach (iii) for a while, and eventually go to 
(i). As a result, |c| can first increase and then decrease as 
we go to long distances. This will be discussed in more 
detail in Sec. VI (sec Fig. 15). 
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We see that the only completely stable fixed point is 
given by (i). As we approach this point at large distances, 
b and c go to zero as 

b, c ~ L- a , (17) 

while the ratio b/c approaches a constant. 

Let us now consider the three systems shown in Figs. 
1-3. In all the systems, there are four length scales of 
interest. First, there is the microscopic length scale d 
which will be assumed to be much smaller than all the 
other length scales. Then there is the length L$ of the 
various sub-systems, such as the stub in Fig. 1, each arm 
of the ring in Fig. 2 and the fifth wire in Fig. 3. Next, 
we have the thermal coherence length Lp defined as 



where T is the temperature. As mentioned before, we will 
be interested in three different regimes, namely, the ratio 
Lt/Ls being much smaller than 1 (high temperature), 
comparable to 1 (intermediate temperature), and much 
larger than 1 (low temperature). Finally, we have the 
length Lw of the long wires, namely, wires 1 and 3 in 
Figs. 1 and 2, and wires 1, 2, 3 and 4 in Fig. 3. We 
will assume that Lw is much longer than both Ls and 
Lp. The long wires will be assumed to be connected to 
some reservoirs beyond the distance Lw- However, we 
will not need to consider the reservoirs explicitly in this 
paper, and the length scale Lw will not enter anywhere 
in our calculations. 

The interpretation of Lt is that it is the distance be- 
yond which the phase of an electron wave packet becomes 
uncorrelated. This can be understood as follows. We 
recall that if the bias which drives the current through 
a QW system is infinitesimal, then the electrons com- 
ing into the QW from the reservoirs have an energy 
Ep = h 2 k F /(2m), where Ep is the Fermi energy in the 
reservoirs. At a temperature T, the electron energy will 
typically be smeared out by an amount of the order of 
ksT. The uncertainty in energy is therefore given by 
ksT = AE = hvpAkp, where we have used Eq. (4). 
Hence, Afc^ = kpjT '/ '(hvp) = 1/Lp. If an electron with 
one particular wave number kp travels a distance L, the 
phase of its wave function changes by the amount kpL. 
Hence, the phases of different electrons whose wave num- 
bers vary by an amount Akp will differ by about tt (and 
can therefore be considered to be uncorrelated) if they 
travel a distance of about ttLt- Hence Lp (or irLp) can 
be thought of as the phase relaxation length of a wave 
packet [30]. 

We can now discuss in broad terms the RG procedure 
that we will use for the various systems. In each case, 
we will begin at the microscopic length scale d with cer- 
tain values for the entries of the 3x3 S-matrices at the 
various junctions. We will use Eq. (6) to evolve all the 
S'-matrices. We will follow this evolution till we get to 



the length scale Ls or Lp, whichever is shorter. Two 
possibilities arise at this stage. 

(i) If Lp is less than Ls, we will stop the RG flow at 
the length scale Lp, and then calculate the transmission 
probabilities as discussed in Sec. III. 

(ii) If Lp is larger than Ls, we will stop the RG flow of the 
3x3 matrices at the length scale Ls- Much beyond that 
length scale, the various systems shown in Figs. 1-3 look 
rather different since it no longer makes sense to consider 
the different junctions (and their S'-matrices) separately. 
In particular, the stub and the ring systems look like two 
long wires joined at one point, while the four- wire system 
looks like four long wires joined at one point. Thus they 
all look like systems with only one junction as indicated 
in Fig. 4. This junction is described by an effective 
S- matrix which is 2 x 2 for the stub and ring systems, 
and 4 x 4 for the four-wire system. As we will discuss 
for the different systems in Sees. IV- VI, the effective S- 
matrix is obtained by appropriately combining the 3x3 
S-matrices at the various junctions at the length scale 
Ls; we can think of this process as "integrating out" the 
sub-systems of length Ls- Then we will continue the 
RG flow beyond the length scale Ls, but now with the 
effective S-matrices. This will continue till we reach the 
length scale Lp. At that point, we stop the RG flow and 
compute the transmission probabilities as shown in Sec. 
III. 

[The reason for stopping the RG flow at Lp in all cases 
is that the amplitudes of the various Friedel oscillations 
and the reflections from them (caused by interactions) 
and from the junctions are not phase coherent with each 
other beyond that length scale. Hence all these reflec- 
tions will no longer contribute coherently to the renor- 
malization of the scattering amplitudes described by the 
various S-matrices.] 
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(a) (b) 

FIG. 4. Effective descriptions of the various systems at low 
temperature, Lp » Ls- The stub and ring systems effec- 
tively reduce to a two- wire system with a junction as in (a), 
while the four-wire system reduces to a four-wire system with 
a junction as in (b). 

To summarize, we will carry out the RG flow in one 
stage from the length scale d up to the length scale Lp, 
if Lp < Ls- If Lp > Ls, we will study the RG flow 
in two stages; the first stage will be with one kind of 
S- matrix from d to Ls, while the second stage will be 
with a different kind of S-matrix from Ls to Lp. The 
two kinds of S-matrices will be connected to each other 
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at the length scale L$ as discussed in Sees. IV- VI. In 
all cases, when we finally stop the RG flow (after one 
stage or two), we will compute the transmission proba- 
bilities. The procedure for doing this will be discussed in 
the next section. [In reality, we expect a smooth cross- 
over from one stage of the RG procedure to the other at 
some length scale which is of the order oi Ls- For the 
sake of computational simplicity, however, we are adopt- 
ing an RG procedure which changes abruptly exactly at 
Ls- We should also note that the RG procedures that 
we will follow are really only valid for Lt << Ls and 
Lt >> Ls- However, we will assume for convenience 
that it is a reasonable approximation to use the same 
procedures all the way up to Lt = Ls-] 

In all the numerical results presented in Sees. IV- VI, 
we will take the interaction parameter a = 0.2 on all the 
wires, and the ratio Ls/d = 10. We will always begin the 
RG flow at the length scale Lt — d, i.e., Lt/Ls = 0.1. 
The values of r' that we will quote in the different figures 
will be the values at Lt/Ls = 0.1. 

III. LANDAUER-BUTTIKER CONDUCTANCE 

In this section, we will discuss how to calculate the 
conductances of the various systems in the three differ- 
ent regimes of temperature. As mentioned already, we as- 
sume that in each of the systems, the long wires are even- 
tually connected to reservoirs through some contacts. For 
a single channel of spinless fermions, there is a resistance 
of e 2 /h at the contacts [30]. [Although the contacts can 
themselves scatter the fermions [16], we will ignore such 
effects here. We are also assuming that the QWs are free 
of impurities. So the only sources of scattering in our 
systems are the junctions.] 

We take the fermions in all the reservoirs to have the 
same Fermi energy Ep, and the net current on all the 
wires to be zero in the absence of any applied voltage on 
the leads. Now suppose that the voltage in reservoir i 
is changed by a small amount Vf, here z = 1,3 for the 
stub and ring systems, and i = 1, • • • , 4 for the four- wire 
system. For \V,\ much smaller than all the other energy 
scales in the problem, such as Ep and fc^T, the net cur- 
rent flowing into wire i (from reservoir i) will satisfy the 
linear relationship [30,31] 

7 ' = J E T V v i ' ( 19 ) 

3 

where the (for i ^ j) define the various transmission 
probabilities, and Tu+1 denotes the reflection probability 
on wire i. The satisfy certain sum rules. Current 
conservation implies that 

E T « = ( 2 °) 



for each value of j. The condition that each of the cur- 
rents must be zero if all the Vi are equal to each other 
implies that 

E T » = ( 21 ) 

3 

for each value of i. This is equivalent to saying that 
changing all the Vi by the same amount does not change 
any of the currents. Thus, if there are TV wires, only 
TV — 1 of the voltages are independent variables as far as 
the currents are concerned. 

We can compute any of the conductances of the system 
if we know the values of all the in Eq. (19). One way 
to define a conductance is as follows [30,31]. We consider 
two of the long wires, say, i and j; we call these the 
current probes, and the currents at these two wires satisfy 
Ii = —Ij. On all the other wires, we impose the voltage 
probe condition I m = 0; this imposes N — 2 conditions on 
the voltages. These conditions imply that there is only 
one independent variable left amongst all the voltages; 
we can take this variable to be Vk — Vj, where k ^ I (in 
general, k, I may or may not be the same as i, j). We can 
now define a conductance of the form 

G ijM = Vk h Vi > (22) 

In Sees. IV and V, we will consider systems which have 
only two long wires labeled 1 and 3 (see Figs. 1 and 2). 
Eqs. (20-21) then imply the relations 

T13 = X31 = — T11 = — T33 , (23) 

and there is only one conductance to consider, namely, 

Gis,3i - X T i3- (24) 

In Sec. VI, we will consider a system with four long wires 
labeled 1, 2, 3 and 4, with complete symmetry between 
wires 1 and 2 on one hand, and between 3 and 4 on the 
other (see Fig. 3). In this case, we have the relations 



T- 










Tn 




T22 = 


T33 — 


T44 


T12 




T34 , 






T13 




T M = 


T23 = 


T24 


Tn 


+ 


T12 + 


2Ti 3 


= 



(25) 



There are many different conductances one can consider 
in this system; the important ones are as follows. We can 
take wires 1 and 2 to be the current probes, and wires 3 
and 4 to be the voltage probes. We then obtain 

e 2 

Gi2,2i = j ( T12 + T13 ) . (26) 

Alternatively, we can take wires 1 and 3 to be the current 
probes, and wires 2 and 4 to be the voltage probes. We 
then find 
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We now have the problem of determining the values of 
T{j in our systems. If Lt > Lg, we will see in Sees. IV- VI 
that one can think of each of the systems as effectively 
having only one junction. In that case, is related to 
the entries of the S'-matrix at that junction as follows, 



1 1{ 

and Tu = \n 



for i ^ j 
- 1 . 



(28) 



(Eqs. (20-21) then follow from the unitarity of the S- 
matrix). On the other hand, if Lt < Lg, then we have 
to consider all the junctions in the system, and the cal- 
culation of Tij involves combining the effects of several 
S-matrices in some way. 

All the three systems of interest to us have two junc- 
tions. [Note that the stub also has two junctions, i.e., 
a three-wire junction at the lower end A, and a one- 
wire junction at the upper end B where we will take the 
S-matrix to be equal to —1 which corresponds to the 
hard wall boundary condition.] An electron which enters 
through one of the long wires has the possibility of bounc- 
ing back and forth many times between the two junctions. 
After several bounces, the electron can emerge from the 
same long wire or from a different long wire. From Figs. 
1-3, we can see that two waves which emerge from the 
system after m and n% bounces will have a difference in 
path lengths which is equal to 2\n\ — n2\Lg. Now we see 
that there are two regimes of temperature which will give 
different answers for the probabilities . If the thermal 
length Lt (which is the phase relaxation length as argued 
earlier) is much smaller than Lg, then the two waves will 
be phase incoherent if ri\ ^ n%. In this case, the con- 
tributions of the two waves to must be added inco- 
herently. On the other hand, if Lt is much larger than 
2 1 rii — n2\Ls 1 then the contributions of the two waves 
will add up coherently. In between these extremes is an 
intermediate regime in which Lt is comparable to Lg] in 
that case, we have only partial coherence, and the two 
waves becomes more and more incoherent as |m — 712I 
increases. 

It is useful to have an expression for the which can 
interpolate all the way from the coherent regime (low 
temperature) to the incoherent regime (high tempera- 
ture). To obtain such an interpolating formula, we use 
the idea of partial coherence caused by phase random- 
ization by a voltage probe which was introduced in Refs. 
[28,29]. We will first summarize this idea, and then de- 
scribe how it can be extended to our problem. 

Consider Fig. 5 in which there is a wire with two ends 
labeled as 1 and 2. At some point labeled A in the middle 
of the wire, there are two wires P\ and P2 which are 
voltage probes, i.e., the net outgoing currents at each 



A 



P 

2 



FIG. 5. A wire with two voltage probes Pi and P2 at a 
point labeled A. The two probes cause phase randomization 
of right and left moving waves respectively. 



of these wires is zero. The four-wire junction at A is 
governed by an S'-matrix of the form 







V^p 
VP 
Vp 



p 





-VP 



-Vp 


VT^P 







(29) 



VT^PJ 



where the columns and rows carry the indices 1, 2, Pi 
and P2 in that order, and p is a real parameter which lies 
in the range < p < 1. (If p = 0, the voltage probes have 
no effect, and phase randomization does not occur). A 
wave traveling right from end 1 can go partially out into 
Pi and partially on to end 2; the part which goes out into 
Pi can re-enter the wire and go on to 2. We now impose 
the phase randomization condition that the wave which 
goes out into Pi can re-enter the wire with an arbitrary 
phase change of <f>i , but it must have the same magnitude; 
this is necessary to ensure the zero current condition since 
Pi is a voltage probe. Upon solving a problem in which 
there is an incoming wave of unit amplitude at end 1, 
and no incoming wave at end 2, we find that there is 
complete transmission of the wave across the point A. 
To be explicit, we take the incoming and outgoing waves 
at A to be of the forms 



4'i 



( 1 \ 



V J 



and tpo 



( \ 

t 

w 



(30) 



respectively. We then use the scattering matrix in Eq. 
(29) to relate ipi and ipo- On eliminating tp, we find that 
the transmission amplitude across A is given by 



t 



1 



V 1 ~P e " 



1 - V^-~P e ' ; 



(31) 



so that |i| = 1. When we calculate any physical quantity 
(such as a transmission or reflection probability) , we will 
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integrate over <j>i from to 2tt. The following identity 
will prove to be useful, 

r = (i-p) h/2 (32) 

Jo 2 7T 

for any integer n. As shown below, the integration over 
4>i reduces the coherence of a wave moving to the right 
from 1 to 2. Similarly, we can introduce a phase change 
of (f>2 for a wave which leaves and re-enters the wire at 
the probe P 2 ; integrating over <p 2 reduces the coherence 
of a wave moving to the left from 2 to 1 . For both right 
and left moving waves, the degree of coherence depends 
on the value of the parameter p as we will now see. 

We consider two waves which travel n\ and n 2 times 
respectively through the point A from left to right. Let us 
suppose that their amplitudes are a\ and a 2 respectively 
in the absence of phase randomization, i.e., for p = 0. 
In the presence of phase randomization, their amplitudes 
will be ait™ 1 and a 2 i™ 2 respectively, where t is given in 
Eq. (31). If these two waves contribute to a transmission 
probability T, the cross-term coming from their interfer- 
ence will be given by a 1 a 2 t™ 2- " 1 + a 2 ait" 1- ™ 2 . We now 
integrate this expression over the variable <j)\. Using Eq. 
(32), we find that 

C ~W~ [ «i a 2 t n2 - ni + a* 2ttl t ni ~ n2 } 
Jo 2?r 

= ( a\a 2 + a; ai ) (l-p)l«i-" 2 l/ 2 . (33) 

We thus see that the phase randomization has the effect 
of multiplying the interference of two terms by a factor 
which interpolates between 1 (i.e., complete phase coher- 
ence between the two waves) for p — and (i.e., no 
phase coherence) for p = 1. Further, the interpolating 
factor depends exponentially on \n\ — n 2 | which is pro- 
portional to the difference between the path lengths of 
the two waves. 

Now we have to implement this idea in the systems 
of interest to us. We do this by generalizing the idea of 
phase randomization at a single point to phase random- 
ization at a continuum of points. Let us assume that 
the density of such points in a wire is given by 11 /Lt, 
where /U is some dimensionless number (which is inde- 
pendent of any temperature or length scale), and that 
the parameter p is the same at each of those points. This 
assumption for the density is motivated by our identifi- 
cation of Lt as the phase relaxation length; the smaller 
the value of Lt, the more frequently phase relaxation 
should occur as an electron travels along the wire. Thus 
the number of phase relaxation points in an interval of 
length L is equal to fiL/Lx- Following arguments sim- 
ilar to the one described above, one can show that the 
interference of two waves which pass through that length 
interval n\ and n 2 times will get multiplied by the fac- 
tor (1 - p )\ni-n 2 \tiL/{2L T )^ Tf we write (1 _ p y/2 = e ~ v , 

where v is a positive dimensionless number, we see that 



the interference between two waves whose path lengths 
differ by AL — \n\ — n 2 |L will carry a factor of 

F = exp [ - v ^ ] . (34) 

Lt 

The high temperature limit {Lt — > 0) corresponds to the 
incoherent case in which we ignore the interference be- 
tween paths with any finite length difference; namely, we 
simply add up the probabilities contributed by different 
paths. 

In our calculations of the transmission probabilities de- 
scribed in Sees. IV- VI, we will require an interpolating 
factor as in Eq. (34) only for the case L T < L s . It 
is only in that regime that our systems have more than 
one junction which allows for a number of different paths 
between any pair of long wires. For Lt > Ls, each of 
our systems effectively simplify to a system which has 
only one junction and, therefore, only one possible path 
between any pair of long wires. Hence there will be no 
need to consider any interference terms for Lt > Ls- In 
order to make our expressions for the transmission prob- 
abilities match as we approach Lt = Ls from above and 
below, we will use an interpolating factor F which is 1 
at Lt = Ls- We will therefore use a formula which is 
motivated by the expression in Eq. (34) (with v set equal 
to 1), but which is somewhat modified so that it is goes 
to 1 as Lt approaches Ls from below. We will use the 
following prescription, 

F = exp [ ^ - ^ ] for L T < L s , (35) 
Ls Lt 

for the factor multiplying the interference of two paths 
differing in length by AL. 

IV. THE STUB SYSTEM 

We will now use the ideas developed in the previous 
two sections to study the transmission probability of the 
stub system shown in Fig. 1. The 3x3 scattering matrix, 
called Ssd, which governs the junction labeled A will be 
assumed to be of the form given in Eq. (12-13), with 
complete symmetry between the two long wires labeled 
1 and 3. At the other end of the stub labeled B, we 
will assume a hard wall boundary condition, i.e., perfect 
reflection with a phase change of —1. 

We first consider the RG flow of the transmission prob- 
abilities Tij. There is only one independent quantity to 
consider in this system, namely, T13; all the others are 
related to it by Eq. (23). As outlined in Sec. II, we start 
from the length scale d and initially use Eq. (14) to see 
how the various entries of 5 3 £> flow as functions of the 
length. If Lt < Ls, we follow this flow up to the length 
scale Lt and then stop there. At that point, we compute 
T13 as explained below. 
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If Lt > Ls, we first use Eq. (14) to follow the flow up 
to the length scale Ls- At that point, we switch over to a 
2x2 scattering matrix S2D which can be obtained from 
the matrix S3D that we get at that length scale from the 
RG calculation. The entries of S2d(Ls) and Ssd(Ls) 
can be shown to be related as follows, 



t 2 TJ 



(5*20)11 = (620)33 — r ' — T~ , 

1 + rrj 

(5*2_d)i3 = (620)31 = t' — - 

1 + rrj 



where 



1 



J2k F L s 



(36) 



(37) 



Eq. (36) will be derived in the next paragraph. [The 
phase factor r\ appears because the electrons are assumed 
to have a momentum of Uf in all regions; hence the wave 
functions have factors of exp(ikpx).] Having obtained 
S2D at the length scale Ls, we then continue with the 
RG flow of that matrix following Eq. (10). This flow 
is stopped when we reach the length scale Lt- At that 
point, we compute T13 as explained below. 

Eq. (36) can be derived in one of two ways. The first 
way is to assume an incoming wave with unit amplitude 
on wire 1 and no incoming wave on wire 3, and then 
use the scattering matrix Ssd at junction A and the sign 
change at B. The second way, which is more instructive 
for us and is also easier, is to sum over all the paths that 
an electron can take. For instance, if we consider the 
different paths which go through the stub, we see that 
they are characterized by an integer n = 0, 1, • • • which is 
the number of times a path goes up and down the stub. 
The length of a path which goes from a point just to the 
left of A to itself after going up and down the stub n 
times is given by 2nLs- Summing over all such paths 
leads to the expression 



(£2.0)11 =r' — tr\t + tr/rrjt - trjrrjrr/t + ■ ■ 



t 2 rj 



1 + rr\ ' 



(38) 



which is the first equation in Eq. (36). Similarly, we can 
derive the second equation in Eq. (36) by summing over 
all the paths which go from a point just to the left of A 
to a point just to the right of A 

Let us now calculate the transmission probability T13. 
If Lt < Ls, we have to use S^d to compute an expression 
for T13 with an interpolating factor F as in Eq. (35). 
This is easy to do since we have already found the sum 
over all the paths as in Eq. (38). According to the phase 
randomization idea discussed in Sec. Ill, the interference 
between two paths characterized by integers n\ and 712 
must be multiplied by a factor F = w here 



cxp 



2(1 - *±) 



(39) 



This follows from the prescription Eq. (35) since the dif- 
ference in path lengths is given by AL = 2|m — n2\Ls- 
On summing up all the terms with the appropriate fac- 
tors of /, we find that 



T13 = t' 2 + - 



t z 



- t z 



Vf 



+ 



+ rrjf l + rr]*f 



(40) 



where we have used some of the relations in Eq. (13). 
Eq. (40) is the desired interpolating expression for T12. 
If we set / = (as we must do for Lt « Ls), we get 
the incoherent expression 



T13 = t' 2 + L - = 1 + r> , 



(41) 



which is independent of rj. On the other hand, if we set 
/ = 1 (as we must do at Lt = Ls), we get the coherent 
expression 



T13 = \t' - 



t 2 T] 



1 + rrj 



(42) 



which is just the square of the modulus of (620)13 given 
in Eq. (36). Eq. (40) interpolates between the coherent 
and incoherent expressions depending on the value of /. 

There is a way of directly obtaining the incoherent ex- 
pression in Eq. (41) without summing over paths. We 
will present this derivation here; as discussed in the next 
section, a similar derivation will work for the ring system 
where it is difficult to classify the different paths in a con- 
venient way and therefore to sum over them. The idea is 
to add probabilities (intensities) rather than amplitudes. 
Consider a situation with the following kinds of waves: a 
wave of unit intensity which comes into the system from 
wire 1, a wave of intensity i r which is reflected back to 
wire 1, a wave of intensity i t which is transmitted to wire 
3, a wave of intensity i u which travels up along the stub 
2, and a wave of intensity id down along the stub. [Note 
that the last four waves are actually made up of sums 
of several waves obtained after repeated travels up and 
down the stub; however, we will not need to explicitly 
sum over all those paths in this way of doing the calcula- 
tion. The summation over paths will be implicit because 
we are assuming that i r , i t and i u and id denote the to- 
tal intensities of those four kinds of waves.] Now we use 
the matrix S^d at junction A. This gives the following 
relations between these intensities, 



i r = r' 2 + t 2 i d , 
it = t' 2 + t 2 i d , 
i u = t 2 + r 2 i d . 



(43) 



Similarly, the total reflection at the end B implies that 
iu = id- Putting these relations together, we obtain 



it 



= 1 + r' 



(44) 
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which agrees with the expression in Eq. (41). 

When Lt becomes equal to Ls, T^{Ls) is equal to 
|(5 , 2d)i2| 2 where (£20)12 is given in Eq. (36). [Our for- 
malism is designed to ensure that we get the same value 
of Tig at Lt = Ls whether we approach that point from 
the high temperature or the low temperature side.] Us- 
ing the parametrization in Eq. (9) and the RG equations 
in Eq. (10), we see that A at a length scale Lt > Ls is 
related to its value at the length scale Ls as follows, 

KLt) = ( ^ Y KL S ) . (45) 

Then T 13 (L T ) is given by 1/(1 + A 2 (L T )). 

In the coherent regime given by Lt > Ls, we observe 
that T13 is equal to 1 if r\ = — 1 and if r\ = 1; this follows 
on using Eqs. (13) and (42). We will call these resonances 
and anti-resonances respectively; they arise due to inter- 
ference between the different paths. For these two special 
values of 77, T13 remains stuck at 1 and and does not 
flow under RG. For any other value of 77, T13 starts at a 
value which is less than 1; it then flows towards zero till 
the RG evolution stops at the length scale Lt- Note that 
by changing the electron momentum kp (this can be done 
by changing the gate voltage), we can vary the value of 
77 and therefore of the matrix elements in Eq. (36); we 
can therefore, in principle, tune the system to resonance. 
[This is in contrast to a single wire system with an im- 
purity where one can change the matrix elements of S2D 
only by varying the strength of the impurity potential 
which may not be easy to do experimentally.] 
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FIG. 6. T13 for the stub system as a function of Lt/ Ls for 
a = 0.2, L s /d = 10, 77 = e ilr/2 , and different values of r'. The 
four sets of curves are for r = —0.10, —0.33, —0.54 and —0.80 
from top to bottom. 

In Figs. 6-7, we show T 13 as a function of Lt/Ls for 
various values of r 1 and 77, with a = 0.2 on all the three 
wires, and Ls/d — 10. In Fig. 6, we have considered 
four different values of r'. Of these values, the first one is 
greater than the unstable fixed point value of — 1/3 given 
in Eq. (15), the second is equal to —1/3, and the last 
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FIG. 7. T13 as a function of L T /L S for a = 0.2, L s /d = 10, 
r' — —0.54, and different values of 77. The four sets of curves 
are for 77 = -1, e i0 - 87T , e i0S7T and 1 from top to bottom. 

two are less than —1/3. In the incoherent regime, we see 
that T13 increases in the first case, does not change in the 
second case, and decreases in the last two cases. In Fig. 
7, we show T13 as a function of Lt/Ls for four different 
values of 77. In the coherent regime, we see that T13 
remains stuck at 1 and for rj = — 1 and 1 respectively, 
while it decreases for the other two cases. It is clear that 
modifying 77 (by changing the gate voltage) can lead to 
large changes in T13 in the coherent regime. 

V. THE RING SYSTEM 

We now turn to the ring system shown in Fig. 2. We 
will assume that both the junctions A and B are de- 
scribed by the same 3x3 scattering matrix S^d given in 
Eqs. (12-13), with complete symmetry between the two 
arms of the ring labeled 2 and 4. 

The RG evolution of the transmission probabilities 
is studied in the same way as for the stub system. (Once 
again, there is only one independent quantity to consider 
here, namely, T13; the others are related to it by Eq. 
(23)). We start from the length scale d and initially use 
Eq. (14) to see how the various entries of the two matrices 
S3D flow as functions of the length. If Lt < L5, we follow 
this flow up to the length scale Lt and then stop there. 
At that point, we compute T13 as explained below. 

If Lt > Ls, we first use Eq. (14) to follow the flow 
up to the length scale Ls- At that point, we switch over 
to a 2 x 2 scattering matrix S2D which can be obtained 
from the matrix S3D that we get at that length scale from 
the RG calculation. For the ring system, the off-diagonal 
matrix elements of £2o(L,s) are related to the parameter 
r' appearing in S3d(Ls) as follows [8]. 

(£20)13 — (£20)31 

2cos($/2) 77 1 / 2 (1 - 77) (-2r') (1 + r') 
[1 + (l + 2r')77] 2 - 2(1 + r') 2 (1 + cos$) 77 ' 
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(46) 



where 77 = e F s , and $ is a dimcnsionlcss number 
which is related to the magnetic flux <ps enclosed by the 
ring through the expression 



$ = 



e<pB 
he 



(47) 



Eq. (46) will be derived in the Appendix. Having ob- 
tained S2D at the length scale L$, we continue with the 
RG flow following Eq. (10). The flow is stopped when 
we reach the length scale Lt- At that point, we compute 
X13 as explained below. 

As shown in the Appendix, Eq. (46) can be obtained 
by assuming an incoming wave with unit amplitude on 
wire 1 and no incoming wave on wire 3, and then using 
the scattering matrices at junctions A and B [8]. One 
might think of deriving Eq. (46) by summing over all 
paths which go from wire 1 to wire 3, just as we did for the 
stub system. However, it seems very hard to enumerate 
the set of paths for the ring system in a convenient way. 
This is because there are two arms, and a path can go 
into either one of the two arms every time it encounters 
one of the two junctions. 

This difficulty in summing over paths also makes it 
hard to find a simple interpolating formula for the con- 
ductance in the regime Lt < L$. To see this more 
clearly, we first note that in the stub system, two paths 
which have equal lengths must necessarily be identical 
to each other. Any point which phase randomizes waves 
moving in one particular direction will therefore occur 
the same number of times in the two paths. Hence, the 
interference between the two paths will not come with 
any powers of either the phase factor 77 or the phase ran- 
domization factor /. If there are two paths of unequal 
lengths 1n\Ls and 2n 2 Ls in the stub system, then any 
point which phase randomizes waves moving in a par- 
ticular direction point will occur m times in one path 
and ri2 times in the other path. Therefore the interfer- 
ence between the two paths will come with a factor of 
^m-ns j|ni-n 2 | xhus, the power of 77 and the power of / 
are always related to each other in a simple way. (This is 
why a factor of / always accompanies a factor of 77 or 77* 
in Eq. (40)). The situation is quite different in the ring 
system. Here the powers of 77 and / are not necessarily 
related to each other in any simple way. For instance, 
consider a path which enters the system through wire 1, 
goes into the arm 2 and leaves through wire 3, and a sec- 
ond path which enters through wire 1, goes into the arm 
4 and leaves through wire 3. These two paths have the 
same length Ls; the interference between the two will 
therefore not carry any powers of 77. However, a phase 
randomization point which lies on one path will not lie 
on the other path. Hence, the phase randomizations will 
not cancel between the two paths, and the interference 
between the two paths will carry a factor of /. Thus, 
there is no general relation between the power of 77 and 



the power of /. This makes is difficult to find an inter- 
polating expression for T13. 

Even though we do not have an interpolating expres- 
sion for the ring system, we can obtain an incoherent 
expression for T13 by following a procedure similar to the 
one we used for the stub system. The idea again is to 
add probabilities rather than amplitudes. Consider a sit- 
uation with the following kinds of waves: a wave of unit 
intensity which comes into the system from wire 1, a wave 
of intensity i r which is reflected back to wire 1, a wave of 
intensity it which is transmitted to wire 3, waves of in- 
tensity i 2r and i% which travel respectively from junction 
A to junction B and vice versa along wire 2, and waves of 
intensity i 4r and in which travel respectively from junc- 
tion A to junction B and vice versa along wire 4. [Note 
as before that the last six waves are actually made up 
of sums of several waves obtained after repeated bounces 
from the two junctions. We do not need to explicitly 
sum over all these paths because we are assuming that 
i r ,it,--- denote the total intensities of these six kinds 
of waves.] Now we use the matrices at junctions A and 
B. This gives the following relations between the various 
intensities, 

i r = r 2 + t 2 ( i 2 i + in ) , 
it = t 2 ( i 2r + Ur ) , 



t 2 + r' 2 i 2l + t' 2 ni 
t 2 + r' 2 ni + t' 2 i 2l 



l2r 
l 4r 

12 ■ ,12 ■ 

121 = r l2r + t l ir , 

in = r' 2 i 4r + t' 2 i 2r ■ 



(48) 



Solving these equations and using some of the relations 
in Eq. (13), we find the incoherent expression for T13 
(valid for Lt « L$) to be 



H = - 



2r' (1 + r') 
1 + ? + r' 2 



(49) 



which is independent of both 77 and $. 

At the point L T = L s , T 13 (L S ) is given by |(5 2 d)i3| 2 
where (S 2 d)i3 is given in Eq. (46). We can again use 
the parametrization in Eq. (9) and the RG equation in 
Eq. (10) to we obtain T 13 {L T ) = 1/(1 + A 2 (L T )), where 
X(Lt) is given in Eq. (45). 

In Fig. 8, we show T13 as a function of Lt/Ls for 
various values 77, with a = 0.2 on all the wires, and 
L s /d = 10. 

In the coherent regime given by Lt > Ls, we can 
find the conditions under which there are resonances and 
anti- resonances in the transmission through the ring, i.e., 
T13 = 1 and respectively. (Some of these conditions 
have been discussed in Ref. [8]). We find that Ti 3 = 1 
for the following values of 77, e l * and r' . 
(i) 77 = e 4 * = 1, and r 1 can take any value. Note that 
for these values of 77 and e J * , there are eigenstates of the 
electron which are confined to the ring. 
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FIG. 8. T13 for the ring system as a function of Lt /Ls for 
a = 0.2, L s /d = 10, r' = -0.18, $ = 0, and different values 
of r\. The four coherent curves are for r\ = 1, e l0 08 ' r ! e »o.i67r 
and e 10 ' 4 ^ from top to bottom. The incoherent curve is the 
same for all ij. 
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2kpL s mod 2jt 

FIG. 10. Ti 3 as a function of 2k F L s for a = 0.2, 
Ls/rf = 10, <E> = 0.57T, and different values of Lt/Ls- The 
four sets of curves are for Lt/Ls = 1, 50, 3000, and 10 from 
top to bottom. Note that there is a transmission zero at 
2k F L s = mod 2tt. 
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FIG. 9. Ti 3 as a function oi2k F L s for a = 0.2, L s /d = 10, 
r' = -0.19, $ = 0, and different values of L T /Ls- The three 
sets of curves are for Lt/Ls = 1,50 and 10 5 from top to 
bottom. 

(ii) i] = e ±l *, and r' = 0. For this relation between 
r) and e 4 *, there arc eigenstates of the electron on the 
ring. Further, r' = implies t = which means that 
these eigenstates cannot escape from the ring to the long 
wires. 

(iii) e** = 1, r' = —1/2, and 7] can take any value. Note 
that r' = —1/2 implies r = which means that a wave 
which is coming in on wire 1 (3) suffers no reflection at 
junction A (B). 

[We note that both the numerator and denominator of 
Eq. (46) vanish under conditions (i) and (ii); hence one 
has to take the limit appropriately to see that T13 = 1.] 
Similarly, we hnd that T13 = for the following values 
of r], e l * and r'. 

(i) r/ = 1, and e 4 * and r' can take any values except 1 
and -1 respectively. 



(ii) r' = 0, and rj and e l * can take any values except 
rj = e ±J *. Note that r' = implies that t = which 
means that there is no transmission between the long 
wires and the ring. 

(iii) r' = —1, and r\ can take any value except 1, while 
e l * can take any value. Once again, r' — — 1 implies that 
t = which means that there is no transmission between 
the long wires and the ring. 

(iv) e J * = —1, and (77, r') can take any values except 
(1,-1) and (-1,0). 

As in the stub system, if T13 begins with the value 1 
or at Lt/Ls = 1, it remains stuck there and does not 
flow under RG as we go to larger length scales. For any 
other starting value of T13, it flows towards zero till the 
RG evolution stops at the length scale Lt- It is interest- 
ing to consider the shape of the resonance line which is 
a plot of T13 versus the momentum Uf (or, equivalcntly, 
rj) at very low temperatures. As discussed in the fol- 
lowing paragraph, one finds that the line shape becomes 
narrower with decreasing temperature, with the width 
at half maximum scaling with temperature as T a . Figs. 
9-10 show this feature qualitatively for the resonances of 
types (i) and (ii) described above. In Fig. 10, we see pairs 
of resonances because T13 has maxima at 2kpL$ equal to 
<& and — <& mod 2-7T. Fig. 11 shows the resonance of type 
(iii). Here Ti 3 is close to 1 for a wide range of kp (or 77) at 
Lt/Ls = 1; this is consistent with the resonance condi- 
tion given in (iii) above. We also observe anti-resonances 
(T13 = 0) in Figs. 10 and 11 at 2k F L s = mod 2tt. 

Let us now discuss the resonance line shape in more de- 
tail. Exactly at a resonance, occuring at, say, kp = kp , 
T13 is equal to 1, and it remains stuck at that value 
no matter how large Lt is. We can now ask: what 
is the shape of the resonance line slightly away from 
kp = kp ? If one deviates from kp$ by a small amount 



12 




2kpL s mod 271 



FIG. 11. T13 as a function of 2k F L s for a = 0.2, 
L s /d = 10, $ = 0, and different values of L T /L S - The 
four sets of curves are for Lt/Ls = 1,150,20000 and 10 6 
from top to bottom. Note that there is a transmission zero 
at 2k F L$ — mod 2%. 



[2k F L s mod 2jc] * (L 1 /L s ) a 

FIG. 12. Ti3 as a function of the scaled variable \2kpLs 
mod 2tt] {Lt/LsT for a = 0.2, L s /d = 10, r' = -0.32, 
$ = 0, and different values of Lt/Ls- The four sets of curves 
are for L T /Ls = 1,2.7,7 .4 and 20 from top to bottom. The 
inset shows the same plots without scaling, i.e., T13 as a func- 
tion of 2kpLs- 



Afc = \kp — fcpo| which is fixed, one finds that the trans- 
mission T13 differs from 1 by an amount of order (Afc) 2 . 
(An example of this is discussed below) . Comparing this 
with the form in Eq. (9), we see that A <~ Afc at the 
length scale L5. Eq. (10) then implies that A will grow 
as Ak(Lp/Ls) a at low temperature; hence T13 will ap- 
proach zero as 1/A 2 <~ T 2a at very low temperature, 
if Afc is held fixed. On the other hand, the width of 
the resonance line at half the maximum possible value 
of T13 is given by the condition that Ak(Lp / L$) a ~ 1, 
which implies that Afc ~ T a . Thus the resonance line 
becomes narrower with decreasing temperature, with a 
width Afc which vanishes as T a . To summarize, T13 de- 
pends on the variables Afc and T through the combina- 
tion x = Ak/T a , and T 13 (x) ~ 1/x 2 as x -> 00. (This 
agrees with the expression given in Rcf. [13] for small 
values of a). As a specific example, let us consider the 
resonance of type (i). We set e J * = 1 and take the limit 
kp — > kp = nn/Ls in Eq. (46). We find that at the 
length scale Ls, 



T13 = 1 - 



(l + 2r') 2 
16r' 2 {l + r'f 



(2AfcL s ) 5 



(50) 



up to order (Afc) 2 . Eq. (10) then implies that at the 
length scale Lp, 



T13 



1 + 



(1 + 2/) 



/\2 



16r' 2 (1 + r') 2 



{2AkL s f{^f a 
Ls 



(51) 



Thus, if r' is held fixed and T13 is plotted against 
Ak(Lp/ Ls) a , we should get the same curve for different 
values of Lt/Ls, provided that the quadratic approxi- 
mation in Eq. (50) holds good. In Fig. 12, we show T13 
as a function of 2AkLs(Lp/ Ls) a for four different values 



of Lp/Ls- We see that the curves agree well with each 
other down to about T13 = 0.7. For comparison, we have 
shown the same plots without scaling in the inset; we see 
that they begin disagreeing below Ti 3 = 0.95. (We find 
similar resonance line shapes in the stub and four-wire 
systems, although we have not shown those plots in this 
paper). 



VI. THE FOUR- WIRE SYSTEM 

Finally, let us consider the four-wire system shown in 
Fig. 3. We will assume that both the junctions A and 
B are described by the same 3x3 scattering matrix S%d 
given in Eqs. (12-13), with complete symmetry between 
the wires 1 and 2 on one side and the wires 3 and 4 on 
the other side. The transmission probabilities enjoy the 
symmetries described in Eq. (25). 

We first consider the RG flow of the transmission prob- 
abilities Tij . Due to the symmetries of the system, and 
the relations in Eqs. (20-21), we see that there arc only 
two independent quantities to consider, namely, Ti 2 and 
T13. Following the formalism in Sec. II, we start from 
the length scale d and initially use Eq. (14) to see how 
the various entries of S 3 d now as functions of the length. 
If Lp < Ls, we follow this flow up to the length scale 
Lp, and then compute T\i and X13. 

If Lp > Ls, we first use Eq. (14) to follow the flow up 
to the length scale L5. At that point, we switch over to a 
4x4 scattering matrix S^d which can be obtained from 
the matrix S%d that we get at that length scale from the 
RG calculation. The entries of S4d(Ls) and S^ptiLs) 
can be shown to be related as follows, 
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(SiD)u =r' + 

(5jo)l2 = t + 
(S , 4L))l3 



irt) 



r 2 r\ 



t 2 rr\ 



1 - 

*v /2 



1 



r 2 r) 



(52) 



/2/,' 



is . (Eq. (52) will be derived in the next 
Having obtained Sad at the length scale Ls, 



where 77 = 1 
paragraph) 

we then continue with the RG flow of that matrix using 
Eq. (6). This flow is stopped when we reach the length 
scale Lt, where we compute T12 and T13. 

As in the stub system, Eq. (52) can be derived in one 
of two ways. The first way is to assume an incoming wave 
with unit amplitude on wire 1 and no incoming waves on 
wires 2, 3 and 4, and then use the scattering matrices 
S3 d at junctions A and B. The second way is to sum 
over all the paths that an electron can take. As in the 
stub system, the different paths going between any two 
of the long wires i and j are characterized by an integer 
n = 0, 1, • • • which is the number of times a path goes 
right and left on the central wire labeled 5. For instance, 
the sum over paths which go from a point on wire 1 lying 
very close to the junction A to itself gives the series 



(S 4 D)n=r' + t^'W^t 



+ trf> 2 rrfl 2 rrf> 2 rrf/ 2 t 



+ 



(53) 



which agrees with the first equation in Eq. (52). Simi- 
larly, we can derive the other expressions in Eq. (52). 

Let us now calculate the transmission probabilities Tyi 
and T13. If Lt < Ls, we have to use Ssd to compute 
expressions for with an interpolating factor / as in 
Eq. (39). This is as easy to do here as in the stub system 
since we know how to explicitly sum over all the paths. 
The interference between the contributions of two paths 
characterized by integers n\ and ti 2 must be multiplied 
by a factor /K""^ where / is given in Eq. (39). On 
summing up all the terms with the appropriate factors of 
/, we find that 



Tia = t' 2 + 



t 2 r 2 



2(1 + r 2 ) 



t 2 (t'r 



vf 



T13 



t 2 



2(l + r 2 )J\l-r 2 T]f \-r 2 rff)' 
t 2 r 2 ( nf | V*f \ 
1 — r 2 ri* f J 



2(1 + r 2 ) 2(1 + r 2 ) Vl-r 2 ??/ 



(54) 



These are the desired interpolating expressions for Ti 2 
and T13. If we set / = (as we must do for Lt « Ls), 
we get the incoherent expressions 



T 12 = t' 2 + 



1J2 



t 2 r 



2(1 + r 2 ) 



■13 



(l + r')(2 + 5r' + 4r /2 ) 

2(1 + 2r> + 2r' 2 ) 
r'(l + r') 



2(1 + r 2 ) 



2(1 + 2r' + 2r> 2 ) 1 



(55) 



which arc independent of 77. On the other hand, if we set 
/ = 1 (as we must do at Lt = Ls), we get the coherent 
expressions which are given by the square of the modulus 
of the entries (£40)12 and (54.0)13 m Eq. (52). 

As in the stub system, there is a way of directly ob- 
taining the incoherent expression in Eq. (55) without 
summing over paths, by adding probabilities rather than 
amplitudes. We consider a situation with the following 
kinds of waves: a wave of unit intensity which comes into 
the system from wire 1, waves of intensity 12, 23 and 14 
which go into wires 2, 3 and 4, and waves of intensity 
i r and %i which travel right and left respectively on wire 
5. We then use the matrices 530 at junctions A and B 
to write down the relations between all these intensities, 
and then solve for i 2 and z 3 . This reproduces the results 
inEq. (55). 




— T 13 Incoherent 
.-. T 13 Interpolating 
-- T Coherent 
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FIG. 13. T13 for the four-wire system as a function of 
Lt j Ls for a = 0.2, L s /d = 10, and 77 = 1. The 
four sets of curves (incoherent and interpolating) are for 
r' = -0.43, -0.28, -0.17 and -0.07 from top to bottom. The 
coherent curve remains stuck at 1/4 for all values of r' since 
77=1. 

If Lt > Ls, T12 and T13 are equal to | (540)12 1 2 and 
|(54o)i3| 2 , where (540)12 and (Sad) 13 are given in Eq. 
(52). In this regime, the RG flow has to be carried out 
numerically for the reasons explained after Eq. (16). In 
general, we find that at long distances, b and c flow to 
zero as indicated in Eq. (17); hence Ti 2 and T13 go to 
zero as L~ 2a . 

In the coherent regime given by Lt > Ls, we observe 
that T12 and T13 are both equal to 1/4 if cither 77 = 
1 or r' = —1/2. We may call these resonances since 
the maximum possible value of T13 which is allowed by 
the form of the matrix in Eq. (16) is 1/4. If 77 = 1, 
T13 remains stuck at 1/4 and does not flow under RG. 
This can be seen in Fig. 13 where we show T13 as a 
function of Lt/Ls for various values of r'. For any other 
value of 77, T13 flows till the RG evolution stops at the 
length scale Lt- (As discussed below, T 13 can sometimes 
increase before eventually decreasing towards zero at very 
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— T Incoherent 






■ -■ T [3 Interpolating 


0.05 




-- T [3 Coherent 



ft.l 1 10 

FIG. 14. Ti3 as a function of Lt/Ls for a — 0.2, 
Ls/d — 10, and r = —0.28. The four sets of curves (inter- 
polating and coherent) are for rj = 1, e i0 - 3lT , e i0 6n and e i0 - 9n 
from top to bottom. The incoherent curve is independent of 
V- 



0.22 1 

0.21- 

■ 
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H= ■ / \ ■ 

0.19- / \ - 

is I — — — — — — 

1 10 100 1000 le+04 le+05 le+06 

V L s 

FIG. 15. Ti3 as a function of Lt/Ls for a = 0.2, 
L s /d = 10, r' = -0.33 and r? = e* - 52 ^. T 13 first increases 
and then decreases at very low temperatures. 

low temperatures). As in the stub system, we can vary 
the value of r? and therefore tune the system to resonance 
by changing the electron momentum hp. 

In Fig. 14, wc show T13 as a function of Lt/Ls for 
various values of 77. In Fig. 15, we show the cross-over 
behavior of T13 mentioned in Sec. II. In the coherent 
regime, for certain ranges of values of r' and 77, T13 first 
increases and then decreases at very low temperatures. 

VII. DISCUSSION 

In this work, we have derived the RG equations and the 
transmission probabilities (and conductances) for three 
systems of experimental interest. The RG flows and the 
consequent power-laws in the temperature and length de- 
pendences of the conductances are purely a result of the 
interactions in the wires; there is no RG flow if the in- 



teraction parameters on are all zero. A peculiarity of our 
RG formalism is that it has two stages which work in the 
regimes of high and low temperature respectively. We 
abruptly switch between the two stages when we cross 
the point Lt/Ls = 1. It would be useful to develop an 
interpolating formalism for the RG flow which can vary 
smoothly across the intermediate range of temperature. 
In our way of deriving the RG equations, this may re- 
quire an analysis of the way in which Friedel oscillations 
from two junctions interfere with each other. 

Our results should be applicable to the systems men- 
tioned earlier such as multi-arm quantum wires [1,2], var- 
ious kinds of carbon nanotubes [3,4], and systems with 
other kinds of geometry [5-7]. While some of the early 
experiments focussed on electronic transport in the pres- 
ence of an external magnetic field and the effects of ge- 
ometry, measuring the various conductances at different 
temperatures (and, if possible, different wire lengths) 
should reveal the interaction induced power-laws dis- 
cussed in our work. Note that a spread in the phase (as 
discussed after Eq. (18) in Sec. II) and phase randomiza- 
tion (as discussed in Sec. Ill) are the only effects of ther- 
mal fluctuations that we have considered in this work. 
We have ignored other effects of finite temperature, such 
as momentum relaxation by inelastic scattering, and cor- 
rections to the Landauer-Buttiker conductances due to 
thermal broadening of the Fcrmi-Dirac distribution near 
the Fermi energy. An application of our work to exper- 
iments would require one to disentangle these other fea- 
tures before the effects of interactions can become visible. 
Typically, the temperature at which these experiments 
are done is about 0.1 — 1° K, while the Fermi energy is 
about 10° K which is much larger [6]; hence the thermal 
broadening effect is expected to be small. 

One limitation of our work is that we have assumed lin- 
ear relations between the incoming and outgoing fermion 
fields. In principle, other interesting things can happen 
at a junction, particularly if we consider the case of spin- 
ful fermions and if some of the wires are superconducting 
rather than metallic. For instance, there may be Andrccv 
reflection in which an electron striking the junction from 
one wire is reflected back as a hole while two electrons 
are transmitted into some of the other wires [5,19]. It 
would be interesting to study these phenomena using the 
techniques developed in this paper. 
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APPENDIX 

We will derive Eq. (46) here. We consider a situation 
with the following kinds of waves: an incoming wave on 
wire 1 whose amplitude is unity just to the left of junction 
A, an outgoing wave on wire 1 whose amplitude is -01, out 
just to the left of junction A, an outgoing wave on wire 3 
whose amplitude is V"3,out just to the right of junction B, 
waves on wires 2 and 4 which are incoming at junction A 
and have amplitudes an d i'tia J us t t° the right of A, 
waves on wires 2 and 4 which are outgoing at junction A 
and have amplitudes V^out an< ^ i'taat J us ^ *° the right of 

A, waves on wires 2 and 4 which are incoming at junction 
B and have amplitudes V'ifin an d V^in J us t to the left of 

B, and waves on wires 2 and 4 which are outgoing at 
junction B and have amplitudes tp B out an d 4>4out j ust to 
the left of B. Our aim is to find an expression for the 
transmitted amplitude 03, out- 

The Schrodinger equation relates many of the ampli- 
tudes introduced above to each other. This is because of 
the following features: a wave which travels a distance x 
picks up a phase of e lkpX (we are assuming that all the 
particles have momentum hp), a wave which travels an- 
ticlockwise around the ring from junction A to junction 
B or vice versa picks up a phase of e**' 2 , and a wave 
which travels clockwise around the ring from junction A 
to junction B or vice versa picks up a phase of e~**/ 2 . 
This gives us the following relations: 

V%in — e %,out ! 

I A _ ikpLs+i^/2 , B 

^in = e ifcFLs - i * /2 < out . (56) 

Now we use the form of the scattering matrices in Eq. 
(12) at the two junctions. At junction A, we have 

V'l.out = r + t ( V^in + i>t,in ) . 
V^out = t + r'^4 M + t'^in , 

< out - t + r'V^ in + tty^in • (57) 

At junction B, we have 

03,out = t ( V^in + in ) . 
^out = r V^in + t'^in 7 

V'fout = r V£ in + *V^in • (58) 

Using Eqs. (56-58), we obtain the expression for 
(£20)31 = 03, out given in Eq. (46). 
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